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Effect  of  gas  injection  on  transition  in 
hypervelocity  boundary  layers 


J.S.  Jewell1, 1.A.  Leyva2,  N.J.  Parziale1,  and  J.E.  Shepherd1 


1  Introduction 

A  novel  method  to  delay  transition  in  hypervelocity  flows  in  air  over  slender  bod¬ 
ies  by  injecting  CO2  into  the  boundary  layer  is  presented.  The  dominant  transition 
mechanism  in  hypersonic  flow  is  the  inviscid  second  (Mack)  mode,  which  is  associ¬ 
ated  with  acoustic  disturbances  which  are  trapped  and  amplified  inside  the  boundary 
layer  [8].  In  dissociated  CCE-rich  flows,  nonequilibrium  molecular  vibration  damps 
the  acoustic  instability,  and  for  the  high-temperature,  high-pressure  conditions  as¬ 
sociated  with  hypervelocity  flows,  the  effect  is  most  pronounced  in  the  frequency 
bands  amplified  by  the  second  mode  [3].  Experimental  data  were  obtained  in  Cal¬ 
tech’s  T5  reflected  shock  tunnel.  The  experimental  model  was  a  5  degree  half-angle 
sharp  cone  instrumented  with  80  thermocouples,  providing  heat  transfer  measure¬ 
ments  from  which  transition  locations  were  from  turbulent  intermittency  based  upon 
laminar  and  turbulent  heat  flux  correlations.  An  appropriate  injector  was  designed 
and  fabricated,  and  the  efficacy  of  injecting  CO2  in  delaying  transition  was  gauged 
at  various  mass  flow  rates,  and  compared  with  both  no  injection  and  chemically 
inert  Argon  injection  cases.  Argon  was  chosen  for  its  similar  density  to  CO2.  At 
an  enthalpy  of  approximately  10  MJ/kg  (Eckert’s  reference  temperature  T*  =  2550 
K),  transition  delays  in  terms  of  Reynolds  number  were  documented.  For  Argon 
injection  cases  at  similar  mass  flow  rates,  transition  is  promoted. 


2  Acoustic  delay 

Turbulent  heat  transfer  rates  can  be  an  order  of  magnitude  higher  than  laminar  rates 
at  hypersonic  Mach  numbers.  A  reduction  in  heating  loads  by  keeping  the  bound- 
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ary  layer  laminar  longer  means  less  thermal  protection  is  necessary,  and  hence  less 
weight  to  carry,  or  conversely  more  payload  deliverable,  for  a  given  design. 

The  theory  of  how  energy  from  acoustic  disturbances  is  absorbed  by  relaxation 
processes  is  treated  in,  among  others,  Vincenti  and  Kruger  [12]  which  provides  the 
simplified  Landau-Teller  theory, 

T  =  Ccxp(K2/T)l/3/p ,  (1) 

where  C  and  K2  are  constants,  which  implies  that  the  vibrational  relaxation  time 
T  decreases  with  increasing  pressure  and  temperature.  This  suggests  that  increased 
temperature,  if  pressure  is  constant,  should  generally  permit  acoustic  absorption  at 
higher  frequencies,  as  shown  computationally  by  Johnson  et  al.  [6]. 

Fujii  and  Hornung  [3]  computed  sound  absorption  spectra  for  flows  of  both  air 
and  COt.  While  the  air  flow’s  sound  absorption  curve  peak  occurs  at  much  lower 
frequencies  than  the  calculated  amplification  peak,  in  CO2  the  broad  sound  absorp¬ 
tion  peak  coincides  with  the  calculated  amplification  peaks.  This  coincidence  is 
most  pronounced  at  enthalpies  of  approximately  10  MJ/kg.  Thus,  for  flows  around 
this  enthalpy,  we  might  expect  increasing  the  fraction  of  COt  in  the  boundary  layer 
to  lead  to  significant  acoustic  damping,  and  therefore  delay  transition. 


3  Experiments 

The  facility  used  in  all  experiments  in  the  current  study  was  the  T5  hypervelocity 
reflected  shock  tunnel  at  the  California  Institute  of  Technology;  see  [5]  and  [4].  The 
model  employed  for  all  experiments  was  a  sharp  slender  cone  similar  to  that  used 
in  a  number  of  previous  experimental  studies  in  T5.  It  is  a  5  degree  half-angle  alu¬ 
minum  cone,  lm  in  length,  and  is  composed  of  three  sections:  a  sharp  tip  fabricated 
of  molybdenum  (to  withstand  the  high  stagnation  heat  fluxes),  a  mid-section  con¬ 
taining  a  porous  gas-injector  section  (interchangeable  with  a  smooth,  non-porous 
injector  section  for  control  shots),  and  the  main  body  instrumented  with  a  total  of 
80  thermocouples  evenly  spaced  at  20  lengthwise  locations.  These  thermocouples 
have  a  response  time  on  the  order  of  a  few  /is  and  have  been  successfully  used 
for  boundary  layer  transition  location  in  Leyva  et  al.  [7]  among  many  others.  For  a 
complete  description  of  the  thermocouple  design,  see  Sanderson  [11].  The  conical 
model  geometry  was  chosen  because  of  the  wealth  of  experimental  and  numerical 
data  available  with  which  to  compare  the  results  from  this  program.  A  photograph 
of  the  cone  model  is  shown  in  Figure  1.  The  porous  injector  section  is  4.13  cm  in 
length  and  consists  of  sintered  316L  stainless  steel,  with  an  average  pore  size  of  10 
microns.  A  detail  view  of  the  tip  and  porous  injector  section  is  shown  in  the  bottom 
of  Figure  1.  The  porous  injector  design  was  chosen  in  favor  of  various  injection 
schemes  with  macroscopic  holes,  all  of  which  were  found  to  trip  the  boundary  layer 
and  lead  to  immediate  or  near-immediate  transition.  The  goal  of  the  porous  injector 
was  to  achieve  more  spatially  uniform  injection  flow,  as  discussed  in  Leyva  et  al. 
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[7].  The  injector  section  is  mounted  around  a  plenum  which  supplies  gas  from  a 
tank  instrumented  with  a  pressure  transducer  used  to  compute  mass  flow  rate. 


Fig.  1  Left:  Aluminum  cone,  lm  in  length,  instrumented  with  80  thermocouples  in  20  rows.  Right, 
from  right  to  left:  molybdenum  tip,  plastic  holder  with  316L  stainless  steel  10  micron  porous 
injector,  aluminum  cone  body. 


A  total  of  16  shots  in  T5  make  up  the  data  set  for  the  present  study.  All  were 
nominally  intended  to  have  the  same  flow  conditions:  air  at  10  MJ/kg  and  55  MPa 
in  the  reflected  shock  region.  The  measured  tunnel  conditions  are  presented  in  Table 

1. 


Table  1  Similar  tunnel  conditions  (nominally  10  MJ/kg,  55  MPa),  and  varied  gas  injection  condi¬ 
tions,  with  resulting  transition  Reynolds  numbers,  determined  with  the  intermittency  method. 


H0  [MJ/kg] 

po  [MPa] 

m  [g/s] 

Re, 

2587 

9.94 

51.9 

8.1  (Ar) 

2.09e6 

2589 

10.37 

56.3 

9.3  (C02) 

4.30e6 

2590 

9.92 

55.9 

11.6  (C02) 

4.59e6 

2591 

9.49 

53.2 

4.6  (C02) 

4.23e6 

2592 

9.89 

52.5 

6.9  (C02) 

4.35e6 

2593 

9.99 

55.2 

13.1  (C02) 

4.39e6 

2594 

9.84 

56.1 

16.2  (C02) 

3.54e6 

2596 

10.03 

54.1 

0 

3.88e6 

2597 

10.09 

55.6 

13.9  (Ar) 

3.07e6 

2598 

10.26 

55.3 

0 

3.85e6 

2600 

9.88 

54.7 

11.6  (Ar) 

1.79e6 

2607 

9.80 

54.5 

a 

4.07e6 

2608 

9.75 

55.2 

a 

4.27e6 

2609 

9.92 

55.7 

a 

4.45e6 

2610 

10.37 

54.9 

a 

3.63e6 

2611 

10.44 

54.6 

a 

4.08e6 

a  Solid  plastic  injector  section,  no  flow. 
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4  Results  and  Conclusions 


The  current  study  detected  transition  by  analyzing  the  intermittency  of  the  heat  flux 
signals,  as  described  in  Clark  et  al.  [2]  and  implemented  for  similar  conditions  in 
Mee  and  Goyne  [9].  Intermittency  represents  the  fraction  of  the  run  time  during 
which  flow  over  each  gauge  is  turbulent.  Gauge  signals  are  considered  turbulent 
when  the  signal  is  elevated  above  the  predicted  laminar  value  by  more  than  40%  of 
the  difference  between  the  predicted  laminar  and  turbulent  values.  Gauge  signals  are 
considered  laminar  when  the  signal  is  elevated  above  the  predicted  laminar  value  by 
less  than  20%  of  the  difference  between  the  predicted  laminar  and  turbulent  values. 
For  values  between  20%  and  40%  above  the  laminar  correlation,  the  numerical  in¬ 
termittency  meter  of  Mee  and  Goyne  [9]  is  employed.  As  Mee  and  Goyne  suggest, 
similar  data  sets  (from  the  repeated  shots:  five  smooth  and  two  porous  non-injection) 
are  combined  for  better  intermittency  determination.  Sample  intermittency  plots  for 
two  different  shots  are  presented  in  Figure  2.  Transition  location  is  determined  from 
these  plots  by  noting  where  the  intermittency  trend  departs  zero.  Results  are  pre¬ 
sented  in  the  left-hand  plot  of  Figure  4. 

An  alternate  method  of  determining  transition  location  averages  the  heat  transfer 
rate  from  each  gauge  over  the  entire  test  time  (see  Figure  3).  Transition  is  considered 
to  have  occurred  when  the  trend  departs  the  predicted  laminar  Stanton  number,  as 
in  [1].  Results  for  this  method  are  presented  in  the  right-hand  plot  of  Figure  4. 
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Fig.  2  Top:  Combined  intermittency  by  gauge  location.  Bottom:  The  universal  intermittency  of 
Narasimha  [10],  used  to  find  transition  location. 


Transition  delays  were  documented  in  shots  with  CO2  injection,  compared  both 
to  shots  with  a  porous  injector  but  no  injection,  and  control  shots  with  a  smooth 
injector  section,  as  presented  numerically  in  Table  1  and  graphically  in  Figure  4. 
The  data  show  a  general  trend  of  increasing  delay  with  injection  rate,  before  a  sharp 
dropoff  at  the  highest  injection  rate,  which  may  be  due  to  boundary  layer  detach¬ 
ment.  All  three  Argon  injection  conditions  transitioned  far  earlier  than  any  CO2 
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Heatflux  distribution  -  Shot  #2598 
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Fig.  3  Heat  flux  contour  plots  on  the  developed  cone  surface.  Top  left:  porous  injector  with  no 
injection.  Top  right:  solid  injector  section.  Bottom  left:  Ar  injection  at  1 1 .6  g/s.  Bottom  right:  COi 
injection  at  1 1.6  g/s. 


injection  or  no-injection  conditions.  Average  heat  fluxes  for  several  exemplar  con¬ 
ditions  are  presented  in  Figure  3. 

Future  work  includes  documenting  transition  delay  due  to  gas  injection  at  a 
greater  variety  of  flow  conditions,  as  well  as  more  precise  measurement  of  mass 
flow  rate  and  measurements  of  CCF  mass  fractions  in  the  boundary  layer.  Special 
emphasis  will  be  placed  upon  finding  a  condition  for  which  natural  condition  occurs 
closer  to  the  center  of  the  test  article,  so  that  larger  delays  may  potentially  be  mea¬ 
sured.  Coordination  with  the  numerical  efforts  of  the  G.V.  Candler  group  will  also 
continue,  as  described  in  Wagnild  et  al.  [13]. 
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Fig.  4  Left:  Transition  Reynolds  number,  determined  with  the  intermittency  method,  plotted 
against  injection  mass  flow  rate.  Argon  injection  data,  which  all  lie  far  below  this  Reynolds  num¬ 
ber  range,  are  omitted  for  clarity.  Right:  Transition  Reynolds  number,  determined  with  the  average 
Stanton  number  method,  plotted  against  injection  mass  flow  rate. 


should  not  be  interpreted  as  necessarily  representing  the  official  policies  or  endorse¬ 
ments,  either  expressed  or  implied,  of  AFOSR  or  the  U.S.  Government. 


References 

1 .  Adam,  RH.  (1997)  Enthlapy  Effects  on  Hypervelocity  Boundary  Layers.  PhD  Thesis,  Califor¬ 
nia  Institute  of  Technology,  Pasadena,  CA  91125. 

2.  Clark  J.P.,  Jones  T.V.  and  LaGraff  J.E.  (1994)  On  the  Propagation  of  Naturally-Occurring 
Turbulent  Spots.  Journal  of  Engineering  Mechanics,  Vol.  28,  No.  1:1-19. 

3.  Fujii  K.  and  Hornung  H.G.  (2003)  Experimental  investigation  of  high-enthalpy  effects  on 
attachment-line  boundary  layer  transition.  AIAA  Journal,  Vol.  41,  No.  7. 

4.  Hornung  H.G.  and  Belanger  J.  (1990)  Role  and  techniques  of  ground  testing  simulation  of 
flows  up  to  orbital  speeds.  AIAA  90-1377. 

5.  Hornung  H.G.  (1992)  Performance  data  of  the  new  free-piston  shock  tunnel  at  GALCIT. 
AIAA  92-3943. 

6.  Johnson  H.B.,  Seipp  T.G.  and  Candler  G.V.  (1998)  Numerical  study  of  hypersonic  reacting 
boundary  layer  transition  on  cones.  Physics  of  Fluids,  Vol.  10  No.  10:  2676-2685. 

7.  Leyva  I.A.,  Jewell  J.S.,  Laurence  S.,  Hornung  H.G.  and  Shepherd  J.E.  (2009)  On  the  impact 
of  injection  schemes  on  transition  in  hypersonic  boundary  layers.  AIAA  2009-7204. 

8.  Mack  L.M.  (1984)  Boundary-layer  stability  theory.  Special  Course  on  Stability  and  Transition 
of  Laminar  Flow.  AGARD  Report  709. 

9.  Mee  D.J.  and  Goyne  C.P.  (1996)  Turbulent  spots  in  boundary  layers  in  a  free-piston  shock- 
tunnel  flow.  Shock  Waves,  Vol.  6,  No.  6:337-343. 

10.  Narasimha  R.  (1985)  The  Laminar-Turbulent  Transition  Zone  in  the  Boundary  Layer.  Progress 
in  Aerospace  Sciences,  Vol.  22,  29-80. 

11.  Sanderson  S.  (1995)  Shock  Wave  Interaction  in  Hypervelocity  Flow ,  GALCIT,  California  In¬ 
stitute  of  Technology,  Ph.D.  thesis. 

12.  Vincenti  W.G.  and  Kruger  C.H.  (1965)  Introduction  to  Physical  Gas  Dynamics.  Krieger. 

13.  Wagnild  R.M.,  Candler  G.V.,  Leyva  I.A.,  Jewell  J.S.  and  Hornung  H.G.  (2010)  Carbon  Diox¬ 
ide  Injection  for  Hypervelocity  Boundary  Layer  Stability.  AIAA  2010-1244. 


Distrbution  A:  Approved  for  public  release;  distribution  unlimited 


